Second order quantum phase transition of a homogeneous Bose gas with attractive 

interactions 
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We consider a homogeneous Bose gas of particles with an attractive interaction. Mean field 
theory predicts for this system a spontaneous symmetry breaking at a certain value of the interaction 
strength. We show that at this point a second-order quantum phase transition occurs. We investigate 
the system in the vicinity of the critical point using Bogoliubov theory and a continuous description, 
that allows us to analyze quantum fluctuations in the system even when the Bogoliubov approach 
breaks down. 



I. INTRODUCTION 

Many body problems in a quantum degenerate limit 
can be investigated experimentally and theoretically in 
ultra cold atomic gasses [l|, [2|, 3 ■ ^t the current level of 
the experimental control and detection accuracy various 
phenomena and effects predicted theoretically (for exam- 
ple Mott-Insulator transition [3l) can now be precisely 
investigated and tested. The most fascinating examples 
from our point of view consists of BEC-BCS transition [5| 
and more recently critical behavior at the phase transi- 
tion, including critical exponent of the correlation length 
for a trapped, weakly interacting Bose gas [6|. 

In the present publication we focus our attention on a 
homogeneous Bose gas with attractive interactions in a 
elongated 3D box with periodic boundary conditions. For 
a corresponding ID system the mean field theory (Gross- 
Pitaevskii equation [7| ) predicts a spontaneous symmetry 
breaking when the strength of the attractive interaction 
reaches a critical value Q. On the other hand the one- 
particle density has to be uniform in the ground state of 
the full quantum Hamiltonian, due to the translational 
invariance of the system. The connection between the 
mean field approach, Bogoliubov [9| and the full quantum 
many body theory was analyzed by few authors. For 
instance Ueda and his cooworkers [i3|, considering ID 
Bose gas in the box showed that when the interaction 
strength exceeds the critical value, even weak external 
symmetry breaking perturbation leads to a formation of a 
bright soliton. They also analyzed quantum fluctuations 
(in position and momentum) of the center of mass of the 
soliton [ll| . Finally, in the context of the present study 
we mention the paper of Castin and Hcrzog |12l | where 
they suggested that as long as we restrict the analysis of 
the system to the low order correlation functions it can be 
characterized quite accurately using modified mean-field 
approach. 

Here we consider quantum many body system and we 
show that Hilbert space in the vicinity of the critical 
point can be divided into a subspace, which can be de- 
scribed by the Bogoliubov theory, and a subspace where 



(two) unstable modes (unstable in the Bogoliubov de- 
scription) can be very elegantly treated with the help of 
a so called continuous description. The continuous de- 
scription in our case leads to a Schrodingcr equation of a 
fictitious particle moving in the effective two-dimensional 
potential and a transformation of this potential from 
parabolic to Mexican hat shape reflects the second order 
quantum phase transition. We also use the continuous 
description to analyze fluctuations close to the critical 
point in our homogeneous systems. The analogous anal- 
ysis is also possible for the system in the double well and 
we will present our results (among other things we show 
that in this case there is only one unstable mode) in the 
separate publication [13J . 

The paper is organized as follows. In Sec. [TTl we in- 
troduce the system and shortly present the results of the 
Bogoliubov theory. In Section Hill we restrict the analysis 
to condensate and two unstable modes, using an argu- 
ment (derived in the previous section) that all the other 
modes remain weakly populated. We introduce effective 
Hamiltonian and continuous description, which allows us 
to investigate quantum fluctuations in the system in a 
vicinity of the critical point. We conclude in Sec. lIVI and 
in the Appendix we analyze eigenstates of the system 
and show that for large particle number corrections to 
an effective Hamiltonian are negligible. 



II. BOGOLIUBOV THEORY 

We consider a gas of Bose particles with contact at- 
tractive interaction in a 3D box with periodic boundary 
conditions. The Hamiltonian of the system reads 
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where ip{r) is the bosonic field operator, and g = 
ATrh^a/m, is the coupling constant that characterizes par- 
ticle interaction (a is s-wave scattering length which for 
the attractive interaction is less than zero). In the case 
of homogeneous system it is convenient to switch to a 



momentum basis, 



^(r) 
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where V = L^LyLz is the volume of the box and the sum 



runs over discrete momenta, 

k - 2 f— ^ —\ 
\Lx Ly Lz J 

with integer n^,ny,nz. Then, the Hamiltonian reads 
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Exact diagonaUzation of the Hamiltonian ([4]) is impos- 
sible because dimension of the Hilbert space of the sys- 
tem is enormous. However, we are interested in a ground 
state of the system in the case of weak particle interac- 
tion. Then, one can expect that only one mode macro- 
scopically occupied by atoms and, in the first approxima- 
tion, we may use a mean field approach which relies on 
substitution of the bosonic operator '^(y) by a classical 
field 00 (r). The resulting Gross-Pitaevskii equation 
reveals homogeneous stationary solution (i.e. the con- 
densate wavefunction) 



Uk-Vk = {Uk+Vk) 



(9) 



allows one to write the Hamiltonian ([7]) in a diagonal 
form, Hb = ^k^oEkblik, where the Bogoliubov spec- 



trum reads 
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Within the A^-conserving Bogoliubov theory we can ob- 
tain iV-body ground state in the particle representation 
in a simple form [l|, |l5| 
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and chemical potential fi ~ gN/V. Quantum fluc- 
tuations around the mean field solution can be de- 
scribed within the Bogoliubov theory Q . We employ TV- 
conserving version of the Bogoliubov theory [l4| where 
a part of the Hamiltonian Q, minus a constant term 
liN, with contributions of the second order in Ok (where 
k^O), 



k#0 



2m 
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is substituted by 
Hb = 

where 
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Here 



dNk = {vk\vk). 



(14) 



are eigenvalues of the reduced single particle density ma- 
trix. The total number of atoms depleted from a conden- 
sate mode reads 



dN = Y^ dNk 



(15) 



In the following we will consider the system in a box 
where Lx,Ly < Lz/2. Then, the lowest Bogoliubov exci- 
tation energy corresponds to momenta 



27r , 
±^(0,0,1), 



and it is real provided a parameter 
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where 
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The Bogoliubov transformation 5k = (ufcjAk) 
with 
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is greater than —1. For a = — 1 the energy gap between 
the ground state and the first excited state disappears, 
the homogeneous mean field solution ^ looses its sta- 
bility and the quantum depletion dN diverges because 
dNk^ diverge. More precisely, for a < —1 the mean field 



solution spontaneously breaks the symmetry of the sys- 
tem. In the next section we show that this point can be 
regarded as a critical point. 

We conclude that the Bogoliubov approximation 
breakdowns down at the critical point. The population 
of k-t is diverging, but population of higher momentum 
modes is within this approximation still negligible. In- 
deed, average number of particle in these modes for a 
around —1 can be estimated as follows: 
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and it is not greater than 0.2 in a vicinity of the critical 
point for our box geometry. 

In the next section we analyze higher order terms of 
the Hamiltonian (j4|) and show that one is able to obtain 
a simple and elegant description of the system in a vicin- 
ity of the critical point if the lowest energy modes are 
carefully considered. Strictly speaking we will consider 
only three modes, the condensate mode and k± modes. 



where 



' = (cq + .9^) ( a-|_a+ + a_a^ ) + 5*^ ( S,+a_ + a_^a_ J 
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Notc, that the above Hamiltonian commutes with the 
total momentum operator which in this case is equal to 



27rfi / I 4. 
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Next we represent first excited modes by their symmet- 
ric and antisymmetric combinations defined in Eq. (|12p . 
and describe dynamics of these modes using the position- 
momentum representation 



Xc + ipc , Xs + ips 

Qc = — — a.9 
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III. CONTINUOUS DESCRIPTION 



In this representation the Hamiltonian H' can be split 
into several parts 



In the present section we suggest an effective method 
to describe our system (cold gas of bosonic atoms in a 
box) in a vicinity of the critical point. We consider the 
elongated rectangular box where L^^Ly < Lz/2. From 
the discussion presented above we expect that in this case 
there are two modes corresponding to momenta k± , that 
need careful treatment. 



A. Effective Hamiltonian 

The Hamiltonian restricted to the condensate and k± 
modes reads 

ii « eo I a^a+ + a_a_ I + — I a^a'^a+a- + a_!^a_aoao 
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where we introduced notation dk^ = d±. Since total 
number of particles is equal to N and is conserved we can 
eliminate one of the modes from the above Hamiltonian. 
Here we choose to eliminate the fc = mode. In order to 
do so we substitute operators do and dj by an operator 

N — I d'|_d+ + dLd_ ) (see [l6[). Now our Hamiltonian 
can be written as 



H^^{N'~N)+H', 



(20) 
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The total momentum operator (|2ip is in this representa- 
tion proportional to PcXg — PsXc f^nd it commutes with 
the Hamiltonian Hp (and obviously also with 5H)\ This 
property will be used in the Appendix. 

Up to now we have divided our Hamiltonian into two 
parts and in what follows we will argue that for large 
7V (strictly speaking in the limit of 7V tending to infin- 
ity while a is kept constant) the dominating contribution 
arise from effective Hamiltonian Hp- This is the Hamil- 
tonian of a fictitious particle moving in a two-dimensional 
effective potential. Indeed, upon defining r^ = x\ -\- x'j. 
we can derive a Schrodinger equation 
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where the 2D effective potential is equal to 
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This potential evolves, when we pass through the critical 
point, from a parabolic well (a > —1), through quartic 
well (a = —1) to a Mexican hat shape (a < —1), which is 
a signature of the second order quantum phase transition. 
In the Appendix we present analytic approximations to 
the solution of this Schrodinger equation and estimate 
the contribution from the Hamiltonian 5H in all three 
regions mentioned above. We show that for large N they 
are indeed negligible. 

Our numerical studies are summarized in Fig. [TJ It 
turns out that our estimates from the Appendix are 
in a perfect agreement with the numerical results. To 
prove that the contribution from 5H is small we diago- 
nalize Hamiltonian H' and Hp numerically. In Fig. [1] we 
show energy levels obtained in each case. One can see 
that already for a moderate particle number {N = 500) 
the low-lying energy levels of the two sets coincide and 
for the greater N {N = 10^) even the higher excited 
states are properly reproduced. In conclusion even for 
moderate number of particles Hamiltonian Hp is a very 
good approximation to H' apart from the constant term 
-1- 



3a 

8Af 



B. Critical fluctuations 

Description of a homogeneous Bose gas with attractive 
particle interactions around the critical point becomes 
very simple when wc use the concept of the effective po- 
tential ([24|l . For instance it allows us to study the dynam- 
ical transition of the system when the scattering length 
is changed in time, especially when the system enters the 
critical region. Here, however, we restrict our considera- 
tions to the static (ground state) properties of the system 
around the critical point. Since the effective potential 
undergoes transition from quadratic through quartic to 
Mexican hat, which is a typical signature of the second 
order phase transition, we expect critical fluctuations to 
show up. Hence we search for the observable that will 
show the clear evidence of the critical behavior (will have 
maximal fluctuations at the critical point). 

Let us consider the operator dN representing the num- 
ber of atoms depleted from a condensate. The contribu- 
tion (fT9)) from the modes described by the Bogoliubov 
Hamiltonian is very small and can be neglected, there- 
fore 



dNfunc + fis, (25) 



where 



n, = ^{pl + xl-l), h, = ^{pl+xl-l). (26) 

The depletion, dN = (dN), increases as we approach 
and pass through the critical point, see Fig. [51 In this 
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FIG. 1: (Color online) Energy levels with respect to the 
ground state energy of the Hamiltonian H' (black lines) and 
Hf (red lines). Panel (a) shows results for A'' = 500 while 
panel (b) for iV = 10^ 



figure we compare mean value of atoms depleted from the 
condensate as a function of parameter a for two different 
total number of atoms. Each panel represents two curves 
that practically sit on top of each other; one obtained 
from the Hamiltonian H' and the other from Hp. The 
linear behavior is due to the linear dependence of the 
square of the position of the minimum in the Mexican hat 
— Tq on the interaction strength a, (see the Appendix). 
For large N we can estimate that {fic + fis) — {x^ + 



Now we carry on to the fluctuations. The variance 
of the dN increases as we cross the critical point but if 
we calculate the variance relative to the average value 
it turns out that the critical point region is clearly indi- 
cated by the maximum of such fluctuations. In Fig. [3] we 
present the variable 



AiV: 



{{dN^{dN)Y 
{dN) 



(27) 



as a function of a around the critical point. Wc compare 



again two curves obtained using Hamiltonians H' and Hp 
for two different number of particles. The agreement is 
already satisfactory for 500 particles and it is excellent in 
the case of 10^ particles. We observe that the maximum 
is shifted from the point a = — 1 and this shift tends to 
zero with the increasing number of particles. 

It is interesting to note that for a < — 1 the fluctu- 
ations of dN are much smaller than the fluctuations of 
fic and fis separately. Indeed, for a < —1, the ground 
state solution corresponds to the wavefunction concen- 
trated around a circle with a radius ro [see Eq. ([34])]. and 
because dN is a function of the distance from the origin 
its fluctuations are small. The operators he and hg are 
related to distances from the axes, and their fluctuations 
are larger than the fluctuations of dN if the wavefunc- 
tion is localized around the ring. It means that while 
the number of depleted atoms is expected to be roughly 
the same in each experiment, these atoms may occupy 
differently two orthogonal modes, i.e. the modes p6|) or 
any other orthogonal combination of them. This is the 
origin of spontaneous symmetry breaking and bright soli- 
ton formation, i.e. atoms depleted from a condensate in 
different experiments may differently occupy the modes 
P^ and the resulting density profiles of atomic clouds 
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FIG. 2: (Color online) Mean depletion from the condensate 
as a function of a ior N = 500 (a) and iV = 10^ (b). Black 
continuous lines correspond to the Hamiltonian H' and red 
dashed lines to Hf- 
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FIG. 3: (Color online) Fluctuations of the condensate deple- 
tion (1221) for iV = 500 (a) and N = 10^ (b) as a function of 
a. Black continuous lines correspond to the Hamiltonian H' 
and red dashed lines to Hp- 



IV. CONCLUSIONS 

We consider a ground state of the Bose gas with attrac- 
tive particle interactions in an elongated 3D box. When 
strength of these interactions increases the system passes 
through a critical point, which is indicated by loss of sta- 
bility of the uniform condensate wavefunction and break- 
down of the Bogoliubov approach. In the vicinity of the 
critical point there are two unstable modes and they re- 
quire a careful treatment — all other modes can be ne- 
glected using argument, emerging from Bogoliubov the- 
ory, that they do not get any significant population. 

We introduce an effective hamiltonian and the contin- 
uous description which allows for a simple analysis of the 
unstable modes subsystem. It turns out that the sub- 
system can be effectively described by a 2D Schrddinger 
equation of a fictitious particle moving in an effective 
potential. This potential changes when the interac- 
tion strength increases and it evolves from a parabolic, 
through quartic (at the critical point) to a Mexican hat 
shape which reflects the second order quantum phase 
transition. 

We have analyzed properties of the ground state, en- 
ergy levels structure and quantum fluctuations of the sys- 
tem in a vicinity of the critical point. The fluctuations of 
number of atoms depicted from a condensate (relative to 
an average value) turn out to be maximal at the critical 
point. These fluctuations are smaller than fluctuations of 



numbers of atoms occupying each of the unstable modes 
separately. It means that the depleted atoms may very 
differently occupy two unstable modes in different exper- 
imental realizations. This is an origin of the spontaneous 
symmetry breaking and bright soliton formation that is 
expected to be observed in a single experimental mea- 
surement. 
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V. APPENDIX: EIGENSTATES OF THE 
EFFECTIVE HAMILTONIAN 

Here we present arguments based on the analytic anal- 
ysis to demonstrate that Hamiltonian SH can be ne- 
glected in comparison with Hp for large N when we con- 
sider the ground state of the system. The procedure is 
as follows. First we find a ground state of Hp (we treat 
it separately in the three regions depending on the value 
of a). With some straightforward approximations we are 
able to find analytical expressions for the ground state 
wavefunction and estimate the energy difference AE be- 
tween the ground and the first excited state (with total 
momentum equal to zero, see below). Next we estimate 
the contribution of the Hamiltonian SH using perturba- 
tion theory and show that for large N first order correc- 
tion is much smaller than AE. 

Notice that as we mentioned above total momentum 
operator ([21]) commutes with both Hp and 5H. Hence 
any excitation caused by the latter can not change the 
value of the total momentum, which is equal to zero in 
the ground state. If we use continuous description and 
express the total momentum operator in the polar coordi- 
nates, it turns out that it is proportional to d/d(j). Hence 
any state with total momentum equal to zero do not de- 
pend on (j). In conclusion, to estimate AE we neglect the 
(j) dependence. 

Finally we stress that considerations presented here 
are given only as an example of the back of the envelope 
approximation, since in parallel we perform numerical 
calculations and obtain exact energies and eigenstates of 
both Hamiltonians. 

• The a -|- 1 > case. 

In this case in order to get an analytic expression 
for the ground state we neglect the fourth order 
terms in the effective potential ((24|). Then the 
Schrodinger equation (j23|) can be written in the 



form 



-A^ir) 



1 -\- a 



and the ground state solution is 



Tpir) '^ exp 
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where 



^ho 



l + a 



(28) 



(29) 



(30) 



The difference AE between the energy levels is of 
the order of aj^^ . The fourth order terms omitted in 
Eq. ([^ result in corrections of the order of a\^/N. 
One can omit the fourth order terms in the effective 
potential when af^/N <C a^^ . Inserting the value 
of aiio from the above equation we get [l7[ 



(l + «)^/^»^. 



(31) 



The biggest term in the Hamiltonian SH is of the 
order a\^/N and is smaller than just discussed 
fourth order term. So in the region where the above 
solution is valid the Hamiltonian SH has negligible 
impact. 

The critical point, a = — 1. 

For a = — 1 the Schrodinger equation ([23]) takes 
the form 
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-AMr) + 
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r ip{r) = Eip{r). 
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We switch to the new variable f 
obtain 



r/iVi/6 and 
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(33) 



From dimensional analysis we estimate AE to be 
of the order of N~^/^ . On the other hand the dom- 
inant contribution to the Hamiltonian SH comes 
from the terms Sr' /N and is of order of N~'^/^. So 
in the limit of large N it can be neglected. 

The a + 1 < case. 

For a < —1 the potential in Eq. (|24|) has the shape 
of a Mexican hat. The minimum occurs for 
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N 



8(a-H) 
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(34) 



Around this minimum, we choose the variable R = 
r — To, and rewrite the Schrodinger equation 
as 



EiP{R) 
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i?V(i?), 



use harmonic approximation (sec the comment be- 
low) and finally obtain 

^l^^b{R)-{l + a)R'^P{R)=E^{R). (36) 
The ground state is 

i?2 



ipiR) ~ exp I — 



26L 



where 



bL = 



1 



2(1 + a) 



(37) 



(38) 



Hence, the energy level spacing AE is proportional 
to bj^^ . To justify our approximation leading to the 



harmonic equation (j36p we notice that the omitted 
terms in equation (|35p are of the order r^ b^^ and 



b\^/N . These terms are negligible when 



■ af^ » 



N' 



(39) 



Incidentally it is the same condition as PT|) ob- 
tained earlier. 

Even if the harmonic approximation of equation 
((55|) is valid we still need to estimate the contribu- 
tion of 5H. First order perturbation gives b^^r^/N 
as a leading order. Hence, comparing it with the 
distance between the energy levels, we obtain the 
condition when the Hamiltonian 5H can be ne- 
glected to be 



|l + a| <1 



(40) 
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